List of all input files

/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC1/all_NNNATTAGC_L5Ab_BC1_rnc-HTF_1.unique.hyb.p_values.tab /home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC1/all_NNNCACTAGC_L5Bc_BC1_rnc-HTF_3.unique.hyb.p_values.tab /home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC1/all_NNNGCGCAGC_L5Ac_BC1_rnc-HTF-2.unique.hyb.p_values.tab

list_of_files
$DM_MiSeq_Bc
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/DM_MiSeq/all_NNNCACTAGC_L5Bc_rnc-HTF_1.unique.hyb.p_values.tab"

$DM_MiSeq_Bd
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/DM_MiSeq/all_NNNTCTCTAGC_L5Bd_rnc-HTF_2.unique.hyb.p_values.tab"

$DM_NextSeq_Bc
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/DM_NextSeq/all_NNNCACTAGC_L5Bc_rnc-HTF_1.unique.hyb.p_values.tab"

$DM_NextSeq_Bd
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/DM_NextSeq/all_NNNTCTCTAGC_L5Bd_rnc-HTF_2.unique.hyb.p_values.tab"

$JW_BC1_Ab
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC1/all_NNNATTAGC_L5Ab_BC1_rnc-HTF_1.unique.hyb.p_values.tab"

$JW_BC1_Bc
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC1/all_NNNCACTAGC_L5Bc_BC1_rnc-HTF_3.unique.hyb.p_values.tab"

$JW_BC1_Ac
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC1/all_NNNGCGCAGC_L5Ac_BC1_rnc-HTF-2.unique.hyb.p_values.tab"

$JW_BC4_Bc
[1] "/home/ignatius/PostDoc/2020/StaphCLASH2020/Results/Hybrids_P_Values/JW_BC4/all_NNNCACTAGC_L5Bc_BC4_rnc-HTF_4.unique.hyb.p_values.tab"

list_of_tables <- purrr::map( list_of_files, vroom::vroom )
Rows: 21
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 34
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 26
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 36
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 653
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 202
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 420
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message
Rows: 80
Columns: 37
Delimiter: "\t"
chr [16]: left_hyb_id, left_rna_class, left_chromo, left_name, left_strand, right_hyb_id, right_rna_class, right_chromo, right_name, right_strand, left_transcript_start, l...
dbl [21]: g_x, left_hyb_count, left_genomic_start, left_genomic_stop, g_y, right_hyb_count, right_genomic_start, right_genomic_stop, smaller_max_depth, larger_max_depth, p...

Use `spec()` to retrieve the guessed column specification
Pass a specification to the `col_types` argument to quiet this message

Remove tRNA, rRNA

List of all hybrids

Full join on itself

Overlap same orientation

Given Table of hybrids X and Table of hybrids Y, We join the table with the following rules: * left-hand-side of hybrid in Table X matches with left-hand-side of hybrid in Table Y, and * right-hand-side of hybrid in Table X matches with right-hand-side of hybrid in Table Y

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKCmBgYHtyIGluY2x1ZGUgPSBGQUxTRX0KIyMgTG9hZCBsaWJyYXJpZXMgCmlmKCFyZXF1aXJlKHBhY21hbikpIHsKICBpbnN0YWxsLnBhY2thZ2VzKCJwYWNtYW4iKQogIGxpYnJhcnkocGFjbWFuKQp9CgpwX2xvYWQodGlkeXZlcnNlKQpwX2xvYWQodnJvb20pCnBfbG9hZChzdHJpbmdyKQpwX2xvYWQoaGVyZSkKYGBgCgpgYGB7ciBpbmNsdWRlID0gRkFMU0V9CiMjIERpcmVjdG9yaWVzIE1hbmFnZW1lbnQKYmFzZV9kaXIgPC0gIGhlcmU6OmhlcmUoKSAKZGF0YV9kaXIgPC0gZmlsZS5wYXRoKCBiYXNlX2RpciwgIkRhdGEiKQpyZXN1bHRzX2RpciA8LSBmaWxlLnBhdGgoYmFzZV9kaXIsICJSZXN1bHRzIikKc291cmNlX2RpciA8LSBmaWxlLnBhdGgoYmFzZV9kaXIsICJTb3VyY2UiKQpoeWJfcF92YWxfZGlyIDwtIGZpbGUucGF0aCggcmVzdWx0c19kaXIsICJIeWJyaWRzX1BfVmFsdWVzIikKCmBgYAoKIyMgTGlzdCBvZiBhbGwgaW5wdXQgZmlsZXMgCgovaG9tZS9pZ25hdGl1cy9Qb3N0RG9jLzIwMjAvU3RhcGhDTEFTSDIwMjAvUmVzdWx0cy9IeWJyaWRzX1BfVmFsdWVzL0pXX0JDMS9hbGxfTk5OQVRUQUdDX0w1QWJfQkMxX3JuYy1IVEZfMS51bmlxdWUuaHliLnBfdmFsdWVzLnRhYgovaG9tZS9pZ25hdGl1cy9Qb3N0RG9jLzIwMjAvU3RhcGhDTEFTSDIwMjAvUmVzdWx0cy9IeWJyaWRzX1BfVmFsdWVzL0pXX0JDMS9hbGxfTk5OQ0FDVEFHQ19MNUJjX0JDMV9ybmMtSFRGXzMudW5pcXVlLmh5Yi5wX3ZhbHVlcy50YWIKL2hvbWUvaWduYXRpdXMvUG9zdERvYy8yMDIwL1N0YXBoQ0xBU0gyMDIwL1Jlc3VsdHMvSHlicmlkc19QX1ZhbHVlcy9KV19CQzEvYWxsX05OTkdDR0NBR0NfTDVBY19CQzFfcm5jLUhURi0yLnVuaXF1ZS5oeWIucF92YWx1ZXMudGFiCgpgYGB7cn0KbGlzdF9vZl9maWxlcyA8LSBsaXN0KCAKICBmaWxlLnBhdGgoIGh5Yl9wX3ZhbF9kaXIsICJETV9NaVNlcSIsICJhbGxfTk5OQ0FDVEFHQ19MNUJjX3JuYy1IVEZfMS51bmlxdWUuaHliLnBfdmFsdWVzLnRhYiIgKSwKICBmaWxlLnBhdGgoIGh5Yl9wX3ZhbF9kaXIsICJETV9NaVNlcSIsICJhbGxfTk5OVENUQ1RBR0NfTDVCZF9ybmMtSFRGXzIudW5pcXVlLmh5Yi5wX3ZhbHVlcy50YWIiICksCiAgZmlsZS5wYXRoKCBoeWJfcF92YWxfZGlyLCAiRE1fTmV4dFNlcSIsICJhbGxfTk5OQ0FDVEFHQ19MNUJjX3JuYy1IVEZfMS51bmlxdWUuaHliLnBfdmFsdWVzLnRhYiIgKSwKICBmaWxlLnBhdGgoIGh5Yl9wX3ZhbF9kaXIsICJETV9OZXh0U2VxIiwgImFsbF9OTk5UQ1RDVEFHQ19MNUJkX3JuYy1IVEZfMi51bmlxdWUuaHliLnBfdmFsdWVzLnRhYiIgKSwKICBmaWxlLnBhdGgoIGh5Yl9wX3ZhbF9kaXIsICJKV19CQzEiLCAiYWxsX05OTkFUVEFHQ19MNUFiX0JDMV9ybmMtSFRGXzEudW5pcXVlLmh5Yi5wX3ZhbHVlcy50YWIiICksCiAgZmlsZS5wYXRoKCBoeWJfcF92YWxfZGlyLCAiSldfQkMxIiwgImFsbF9OTk5DQUNUQUdDX0w1QmNfQkMxX3JuYy1IVEZfMy51bmlxdWUuaHliLnBfdmFsdWVzLnRhYiIgKSwKICBmaWxlLnBhdGgoIGh5Yl9wX3ZhbF9kaXIsICJKV19CQzEiLCAiYWxsX05OTkdDR0NBR0NfTDVBY19CQzFfcm5jLUhURi0yLnVuaXF1ZS5oeWIucF92YWx1ZXMudGFiIiApLAogIGZpbGUucGF0aCggaHliX3BfdmFsX2RpciwgIkpXX0JDNCIsICJhbGxfTk5OQ0FDVEFHQ19MNUJjX0JDNF9ybmMtSFRGXzQudW5pcXVlLmh5Yi5wX3ZhbHVlcy50YWIiICkgKSAKCgpuYW1lcyhsaXN0X29mX2ZpbGVzICkgPC0gYyggIkRNX01pU2VxX0JjIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAiRE1fTWlTZXFfQmQiLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICJETV9OZXh0U2VxX0JjIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICJETV9OZXh0U2VxX0JkIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICJKV19CQzFfQWIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIkpXX0JDMV9CYyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSldfQkMxX0FjIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICJKV19CQzRfQmMiICkKCmBgYAoKCmBgYHtyfQoKbGlzdF9vZl90YWJsZXMgPC0gcHVycnI6Om1hcCggbGlzdF9vZl9maWxlcywgdnJvb206OnZyb29tICkKCm5hbWVzKCBsaXN0X29mX3RhYmxlcykgPC0gbmFtZXMoIGxpc3Rfb2ZfZmlsZXMpCgpgYGAKCiMjIFJlbW92ZSB0Uk5BLCByUk5BCmBgYHtyfQoKIyMgTmVlZCB0byBwdXQgdGhlIC4gaW4gYSBwYXJlbnRoZXNpcyAoLikKIyBodHRwczovL3N0YWNrb3ZlcmZsb3cuY29tL3F1ZXN0aW9ucy8zNzc1NzE5MS9waXBlLWEtZGF0YS1mcmFtZS10by1hLWZ1bmN0aW9uLXdob3NlLWFyZ3VtZW50LXBpcGVzLWEtZG90CmZpbHRlcmVkX3RibCA8LSBwdXJycjo6bWFwKCBsaXN0X29mX3RhYmxlcywgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgfnsgKC4pICU+JSBkcGx5cjo6ZmlsdGVyICggbGVmdF9ybmFfY2xhc3MgIT0gInRSTkEiICYKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxlZnRfcm5hX2NsYXNzICE9ICJyUk5BIiAmIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcmlnaHRfcm5hX2NsYXNzICE9ICJ0Uk5BIiAmCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICByaWdodF9ybmFfY2xhc3MgIT0gInJSTkEiICApfSApICU+JQogICAgICAgICAgICAgICAgYmluZF9yb3dzKC5pZD0iRXhwZXJpbWVudF9OYW1lIikgJT4lCiAgICAgICAgICAgICAgICBtdXRhdGUoIGZ1bGxfam9pbl9kdW1teV9pZCA9IDEpCgoKYGBgCgoKIyMgTGlzdCBvZiBhbGwgaHlicmlkcwpgYGB7cn0KCmFsbF9oeWJfaWRzIDwtIGZpbHRlcmVkX3RibCAlPiUKICBkaXN0aW5jdCggRXhwZXJpbWVudF9OYW1lLCBsZWZ0X2h5Yl9pZCApCgpgYGAKCgojIyBGdWxsIGpvaW4gb24gaXRzZWxmCmBgYHtyfQoKZnVsbF9qb2luX3RibCA8LSAgZmlsdGVyZWRfdGJsICU+JQogIGZ1bGxfam9pbiggZmlsdGVyZWRfdGJsLCBieT1jKCJmdWxsX2pvaW5fZHVtbXlfaWQiKSApICU+JQogIGRwbHlyOjpmaWx0ZXIoIEV4cGVyaW1lbnRfTmFtZS54ICE9IEV4cGVyaW1lbnRfTmFtZS55ICkKCmBgYAoKIyMgT3ZlcmxhcCBzYW1lIG9yaWVudGF0aW9uICAKR2l2ZW4gVGFibGUgb2YgaHlicmlkcyBYICBhbmQgVGFibGUgb2YgaHlicmlkcyBZLApXZSBqb2luIHRoZSB0YWJsZSB3aXRoIHRoZSBmb2xsb3dpbmcgcnVsZXM6IAogKiBsZWZ0LWhhbmQtc2lkZSBvZiBoeWJyaWQgaW4gVGFibGUgWCBtYXRjaGVzIHdpdGggbGVmdC1oYW5kLXNpZGUgb2YgaHlicmlkIGluIFRhYmxlIFksIGFuZAogKiByaWdodC1oYW5kLXNpZGUgb2YgaHlicmlkIGluIFRhYmxlIFggbWF0Y2hlcyB3aXRoIHJpZ2h0LWhhbmQtc2lkZSBvZiBoeWJyaWQgaW4gVGFibGUgWQpgYGB7cn0KCnNhbWVfZGlyZWN0aW9uX2NvbXBhcmUgIDwtIGZ1bGxfam9pbl90YmwgJT4lCiAgZHBseXI6OmZpbHRlciggKCBiZXR3ZWVuKCBsZWZ0X2dlbm9taWNfc3RhcnQueSwgIGxlZnRfZ2Vub21pY19zdGFydC54LCAgbGVmdF9nZW5vbWljX3N0b3AueCApIHwKICAgICAgICAgICAgICAgICAgIGJldHdlZW4oIGxlZnRfZ2Vub21pY19zdG9wLnksICAgbGVmdF9nZW5vbWljX3N0YXJ0LngsICBsZWZ0X2dlbm9taWNfc3RvcC54ICkgfAogICAgICAgICAgICAgICAgICAgYmV0d2VlbiggbGVmdF9nZW5vbWljX3N0YXJ0LngsICBsZWZ0X2dlbm9taWNfc3RhcnQueSwgIGxlZnRfZ2Vub21pY19zdG9wLnkgKSB8CiAgICAgICAgICAgICAgICAgICBiZXR3ZWVuKCBsZWZ0X2dlbm9taWNfc3RvcC54LCAgIGxlZnRfZ2Vub21pY19zdGFydC55LCAgbGVmdF9nZW5vbWljX3N0b3AueSApICkgICYKICAgICAgICAgICAgICAgICAoIGJldHdlZW4oIHJpZ2h0X2dlbm9taWNfc3RhcnQueSwgcmlnaHRfZ2Vub21pY19zdGFydC54LCByaWdodF9nZW5vbWljX3N0b3AueCApIHwKICAgICAgICAgICAgICAgICAgIGJldHdlZW4oIHJpZ2h0X2dlbm9taWNfc3RvcC55LCAgcmlnaHRfZ2Vub21pY19zdGFydC54LCByaWdodF9nZW5vbWljX3N0b3AueCApIHwKICAgICAgICAgICAgICAgICAgIGJldHdlZW4oIHJpZ2h0X2dlbm9taWNfc3RhcnQueCwgcmlnaHRfZ2Vub21pY19zdGFydC55LCByaWdodF9nZW5vbWljX3N0b3AueSApIHwKICAgICAgICAgICAgICAgICAgIGJldHdlZW4oIHJpZ2h0X2dlbm9taWNfc3RvcC54LCAgcmlnaHRfZ2Vub21pY19zdGFydC55LCByaWdodF9nZW5vbWljX3N0b3AueSApICkgICkKCgpgYGAKCgoKCgoKCgoKCgoKCgoKCgoK